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^ (57) Abstract: The present invention discloses a method for generating elevation maps or images of a tissue layer/boundary with 
O respect to a fitted reference surface, comprising the steps of finding and segmenting a desired tissue layer/boundary; fitting a smooth 
O reference surface to the segmented tissue layer/boundary; calculating elevations of the same or other tissue layer/boundary relative to 
^ the fitted reference surface; and generating maps of elevation relative to the fitted surface. The elevation can be displayed in various 
O wa y s including three-dimensional surface renderings, topographical contour maps, contour maps, en-face color maps, and en-face 
grayscale maps. The elevation can also be combined and simultaneously displayed with another tissue layer/boundary dependent 
j^- set of image data to provide additional information for diagnostics. 



WO 2007/028531 Al F UK I Illf III! II lllfll lllll IIIN lllff llll f fl f[f Hill IEII! IIIII Iflll (III! 1111 IfilllE III! Illf lllf 



Published: 

with international search report 

For two-letter codes and other ' abbreviations, refer to the "Guid- 
ance Notes on Codes and Abbreviations" appearing at the begin- 
ning of each regular issue of the PCT Gazelle. 



WO 2007/028531 



1 



PCT/EP2006/008465 



METHOD OF BIOIMAGE DATA PROCESSING FOR REVEALING MORE 
MEANINGFUL ANATOMIC FEATURES OF DISEASED TISSUES 

TECHNICAL FIELD OF THE INVENTION 

[0001] One or more embodiments of the present invention relate generally to 

methods for optical imaging of biological samples and for processing such images. In 
particular, the invention is a method for processing a three-dimensional image data set to 
generate elevation maps of tissue layers relative to a fitted smooth surface, which can 
provide more diagnostic information than a pure tissue layer thickness map. Maps of 
elevation may be embodied as three-dimensional surface renderings of elevation, 
topographical maps, or as color or grayscale maps. 

BACKGROUND OF THE INVENTION 

[0002] Measurement of biological tissue surface contour or layer thickness can 

provide useful diagnostic information in various applications. For example, arterial 
plaque thickness is related to the progress of atherosclerosis, carotid vessel wall thickness 
is also an indicator of cardiovascular disease risk; epidermal layer thickness is an 
indicator of burn severity. 

[0003] In ophthalmology, retinal thickness may be abnormally large in cases of 

retinal edema or traction by membranes in the vitreous humor. On the other hand, the 
retina may appear thin in cases of atrophic degeneration, chorioretinitis, or trauma to the 
retina. Meanwhile, changes in retinal thickness may be localized or extend over large 
areas. In certain cases, the overall contour of the retina may become abnormal. For 
example, pronounced myopia, particularly due to posterior staphylomas, may create a 
highly concave retina. Detachment of the retinal pigment epithelium (RPE), subretinal 
cysts, or subretinal tumors may produce a relative convexity of the retina. Therefore, 
mapping the retina contour or retinal thickness makes it possible to determine the extent 
and severity of such conditions and to monitor progress of treatment. 
[0004] In the past, there are a number of well-established biomedical imaging 

techniques that have been used for three-dimensional anatomical mapping of the eye, 
especially the retina, including optical coherence tomography (Zhou, Q. et al. (2004). 
"Mapping retinal thickness and macular edema by high-speed three-dimensional optical 
coherence tomography". Ophthalmic Technologies XIV, SPIE, 5314: 119-125), 
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ultrasound (see for example, US5293871, US5562095), and confocal microscopy (see for 
example, US4838679; R. H. Webb (1996) "Confocal optical microscopy" Rep. Pros. 
EkuL 59 427-471). The three-dimensional data set has also been analyzed to identify 
layered structures in the tissue using a variety of approaches to image segmentation, (see 
5 for example, D. G. Bartsch, et al., (2004) "Optical coherence tomography: interpretation 
artifacts and new algorithm", Proc. SPIE Medical Imaging 2004: Image Processing, 
5370: 2140-2151; H. Ishikawa, et al, (2005) "Macular Segmentation with Optical 
Coherence Tomography". Invest Ophthalmol Vis Sci. ; 46: 2012-201). 
[0005] These prior art methods measured and/or generated a map of a tissue layer 

0 thickness by searching for the borders of the tissue layer structures, figuring out the inner 
and outer boundaries and then finding the distance between the inner and outer 
boundaries. However, a major issue associated with a tissue layer thickness map is that it 
sometimes cannot reveal the diagnostically more meaningful features of a diseased part of 
the tissue. For example, retina thickness is defined as the vertical distance between the 
5 RPE (retinal pigment epithelium) 1 02 and the ILM (inner limiting membrane) 1 04 as 

shown in Fig. 1 . A sharp bump 106 of the retina will often be associated with a rise in the 
RPE 102 as well as the formation of a lesion 108 below the RPE 102, such that the RPE 
also has a broad rise. As a result, a retina thickness map such as the color coded one 
shown in Fig. 2, which corresponds to Fig. 1, cannot reveal the substantially raised bump. 
In fact, the color coded thickness map shows that the thickness will only slightly increase 
near the bump but then return to normal over it. On the other hand, although a 
topographic map or contour of the RPE or ILM may reveal the sharp bump better for this 
illustrated case than the retina thickness map, it would include both the sharp bump and 
the broader warping of the RPE boundary, making it difficult to separate the effect of the 
disease from the overall shape of the RPE boundary 

[0006] In light of the above, there is a need in the art for a method for generating 

elevation maps with respect to a reference fitted surface and for using the reference 
surface as a means of locating three-dimensionally a tissue or a layer or boundary of a 
tissue such as the retina, in order to provide diagnostically more meaningful information 
about potential diseased tissue. 

[0007] The present invention is a novel and non-obvious method wherein a fitted 

reference surface is used to create an elevation map or image of a tissue layer/boundary 
with respect to the fitted reference surface. Use of such a fitted surface can minimize the 



WO 2007/028531 



3 



PCT/EP2006/008465 



perturbations of the surface associated with disease so as to approximate the tissue 
surface that would exist if the tissue were normal. By using such a fitted surface, either 
of the tissue boundary being measured, or a different boundary, the effect of disease or 
injury is isolated from the overall shape of the tissue of interest, providing improved 
5 diagnostic information. In addition to various ways to display the elevation data relative 
to the fitted reference surface, the invention also combines the elevation data with other 
maps or images in order to provide more meaningful information for diagnostics. 

SUMMARY OF THE INVENTION 

1 0 [0008] One or more embodiments of the present invention satisfy one or more of 

the above-identified needs in the art. In particular, one embodiment of the present 
invention is a method for generating elevation maps or images of a tissue layer/boundary 
with respect to the location of a fitted reference surface, comprising the steps of finding 
and segmenting a desired tissue layer/boundary; fitting a smooth reference surface to the 

1 5 segmented tissue layer/boundary; calculating elevations of the same or other tissue 

layer/boundary relative to the fitted reference surface; and generating maps of elevation 
relative to the fitted surface. 

[0009] One aspect of the present invention is to display the elevation in various 

ways including three-dimensional surface renderings, topographical contour maps, 

20 contour maps, en-face color maps, and en-face grayscale maps. 

[0010] Another aspect of the present invention is to combine and hence 

simultaneously display on the same map and/or image two sets of data with one set from 
the elevation relative to a fitted reference surface and the other set from a tissue 
layer/boundary dependent information, including, for example, actual thickness of a 

25 tissue layer, and image signal strength such as reflectance from an OCT system, 

birefringence from a polarization sensitive OCT system or a scanning laser polarimetry 
system, and intensity from a confocal imaging system. 

[001 1] Another aspect of the present invention is to perform the fitting to obtain 

the reference surface in a number of ways, including using a second-order polynomial fit, 
30 or using Zemike or Chebyshev or other polynomials, or Bessel functions, or a portion of a 
sphere or spheroid. Additionally, the fitting can also be performed by excluding certain 
portions of the tissue layer/boundary, i.e. the regions of diseased tissue, from the 
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determination of the fitted reference surface, or fitting on more than one region of the 
tissue layer/boundary or smoothing/filtering a tissue layer/boundary. 
[0012] Still another aspect of the invention is to locate the general tissue 

layer/boundary contour for subsequent scans, which need to follow the tissue contour 
5 closely. 

[0013] Additional aspects of the invention will be set forth in part in the 

description which follows. It is to be understood that both the foregoing general 
description and the following detailed description are exemplary and explanatory only 
and are not restrictive of the invention, as claimed 

10 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0014] Fig. 1 shows OCT images of a retina, illustrating RPE (retinal pigment 

epithelium), the ILM (inner limiting membrane), a sharp bump and a lesion below the 
sharp bump. 

1 5 [0015] Fig. 2 is a color coded retina thickness map corresponding to Fig. 1 

[0016] Fig. 3 shows a flow diagram of the steps of the invented image processing 

method. 

[0017] Fig. 4 shows a retina image with the hue mapped as the square of the 

distance from the fitted RPE reference surface; 
20 [0018] Fig. 5 shows a contour plot using the color coding that represents the 

distance from the ILM to a paraboloid fitted RPE reference surface; 
[0019] Fig. 6 shows a pseudocolor image representing the distance from the ILM 

to a paraboloid fitted RPE reference surface; 

[0020] Fig. 7 shows a three-dimensional rendering of the ILM surface elevation 

25 relative to the paraboloid fitted RPE reference surface, with the color indicating in 
duplicate the same elevation information. 

[0021] Fig. 8 shows a three-dimensional rendering of the actual retina thickness 

superimposed with a pseudocolor image indicating ILM elevation with respect to the 
paraboloid fitted RPE reference surface 
30 [0022] Fig. 9 shows a three-dimensional rendering of the ILM surface elevation 

relative to the paraboloid fitted RPE reference surface superimposed with a pseudocolor 
image indicating actual retinal thickness; 
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[0023] Fig. 1 0 shows a pseudocolor image representing the distance from the RPE 

to a paraboloid fitted RPE reference surface. 

[0024J Fig. 1 1 is a schematic diagram of a basic OCT system capable of 

generating 3D image data that can be used in the method of the subject invention. 

DETAILED DESCRIPTION 

[0025] Fig. 3 shows one preferred embodiment the presently invented method. 

This method is intended to be used on image data obtained from a sample. The 
illustrations in this application are based on image data derived from an optical coherence 
tomography system (OCT) which includes both time domain and spectral domain OCT 
systems. Such an instrument generates 3D intensity data corresponding to an axial 
reflection distribution arising from reflecting features in the eye. As noted above, this 
information is currently used by doctors to view and diagnosis various pathologies in the 
eye. A basic OCT system will be discussed below. 

[0026] Although the illustrated embodiments are limited to OCT data, the image 

processing concepts described herein may be used with 3D image data derived from other 
modalities. For example, image data may be created using various forms of confocal 
microscopy and even ultrasound imaging modalities. 

[0027] The first step in the subject method requires identification of a subset of 

the image data which corresponds to a boundary surface within the sample (step 302). As 
used herein, boundary surface can be a limiting surface in the sample, a surface of a layer 
or other interface. The boundary should have a sufficient linear or 2D extent that it can 
be reasonably fitted to a geometric line or surface. 

[0028] Identification of the image data corresponding to a boundary surface is 

performed using a segmenting function. Methods for finding and segmenting a desired 
tissue layer or boundary surface are well-known in the art. (see for example, H. Ishikawa, 
et al., (2005) "Macular Segmentation with Optical Coherence Tomography". Invest 
Ophthalmol Vis ScL; 46: 2012-201). 

[0029] Once the image data corresponding to the selected surface has been 

segmented, the boundary is fitted to a substantially smooth reference surface (step 304). 
There are a number of well-known methods for fitting a measured surface data points to a 
geometric surface. One example is a second-order polynomial fit. Other functions, 
including Zernike or Chebyshev polynomials, Bessel functions, or a portion of a sphere or 
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spheroid, can also be used for surface fitting. A smooth reference surface can be formed 
by fitting the shape of a tissue layer/boundary with a function of two variables. This 
requires a reasonably accurate segmentation of the chosen tissue layer/boundary and can 
be accomplished using, for example, a low-order polynomial fit in x and y. 
5 [0030] The fitting may encompass the entire tissue layer/boundary or may be 

performed on various regions of the surface, e.g., fitting the foveal pit separately from the 
macula, or excluding pathological regions from the fitting. The reference surfaces can be 
used to define layers in the data that have the retinal tilt and curvature removed. In one 
aspect of the invention, these data points can be used to form en-face images representing 
10 retinal structures in those layers. This presents an advantage over the flat C-scan 

presentation of the data when imaging the curved layers in the anatomy of the eye, since a 
C-scan will only show tangential slices of the retinal layers. 

[0031] As noted above, use of such a fitted surface can minimize the perturbations 

of the surface associated with disease so as to approximate the tissue surface that would 

1 5 exist if the tissue were normal. In this case, the fitting algorithm will function to reject 
points that are associated with the selected boundary surface but exist as a result of the 
disease. By using such a fitted surface, either of the tissue boundary being measured, or a 
different boundary, the effect of disease or injury is isolated from the overall shape of the 
tissue of interest, providing improved diagnostic information. 

20 [0032] In the next step of the method (306), the distances or elevations between 

points on the reference surface and some other feature of interest are calculated. The 
feature of interest may be the actual boundary initially selected so that elevations will 
correspond to the deviations between the selected surface and the associated reference 
surface. The feature of interest can also be another interface or boundary within the 

25 sample. 

[0033] In the next step of the subject invention (308), 2D data sets are generated . 

based on the calculated distances between the reference surface and other feature of 
interest. The 2D data sets can be used to create elevation maps. The elevation maps may 
be created based on pseudocolor or gray scales with elevation encoded as color as shown 
30 in Fig. 6 or intensity (not shown); or as topographical contour maps with elevation 
encoded as contour height as shown in Fig. 5 or as three dimensionally rendered 
topographical maps. These types of maps may also be combined to simultaneously 
display on the same map two sets of data, one for elevation relative to a fitted reference 
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surface and the other for either another elevation relative to another reference surface, or 
an actual tissue layer thickness, or the originally collected image signal strength such as 
OCT or confocal optical signal strength for a tissue layer or other processed/unprocessed 
tissue layer/boundary data such as birefringence measured from a polarization sensitive 
OCT system or a scanning laser polarimetry system. For example, the distance from ILM 
relative to the RPE (i.e. the actual retina thickness) could be displayed as a contour map 
superimposed on a pseudocolor map of ILM elevation relative to a fitted RPE reference 
surface (not shown). Similarly, a pseudocolor map may be applied to a three-dimensional 
surface rendering in order to simultaneously display multiple information on elevation, 
thickness, reflectance or others. 

[0034] In addition to the fact that by analyzing the curvature of the fitted 

reference surface, abnormal tissue layer curvatures (for example, the retina curvature for 
the case of pathologic myopia) can be diagnosed, the present invention has a number of 
other advantages over prior art methods as it can provide additional useful information for 
diagnosing diseased tissues. For example, the fitted reference surface can be used as a 
basis for elevation maps of retinal layers, to diagnose abnormal curvature of the retina, or 
as a guide for subsequent contour-following scans of that eye. Using such a fitted 
reference surface as a basis for "thickness" measurements could give more robust results 
because the exact topography of a deteriorating RPE may be more difficult to determine 
than the general shape of that layer. The reference surface determined by fitting will be 
more consistent than that determined by following the RPE in detail, especially in 
diseased eyes that may have breaks or complex variations in the RPE. As an example, 
map(s) of elevation can be displayed in the form of topographical contour map(s) applied 
to surface renderings of elevation. Fig. 5 shows a color contour coding that represents the 
distance from the ILM to a paraboloid reference surface fitted to the RPE. The corrected 
effective retinal "thickness" relative to the fitted surface is shown in microns on the color 
bar to the right of the map. 

[0035] Additionally, presentation of topographic information relative to a fitted 

reference surface or surfaces can generate images with added information, for example: 
(1) a 2-D false color image giving an en-face presentation of distance from the ILM to 
a reference surface fitted to the RPE can provide information on the effective 
retinal thickness which does not include thickness variations caused by small 



WO 2007/028531 



8 



PCT/EP2006/008465 



perturbations in the RPE. Fig. 6 shows a pseudocolor image representing the 
distance from the ILM to a paraboloid fitted RPE reference surface; 

(2) a 2-D false color image giving an en-face presentation of distance from the actual 
RPE to a reference surface fitted to the RPE itself can highlight localized 

5 variations in the RPE which may be associated with disease. Fig. 10 shows a 

pseudocolor image representing the distance from the RPE to a paraboloid fitted 
RPE reference surface. 

(3) a color mapping for 2-D images or translucent 3-D renderings, in which 
brightness represents reflectance and hue represents distance from the reference 

1 0 surface, can provide an illustration of the height of brightly reflecting layers from 

the reference surface without detailed segmentation of all the layers. For cases of 
complicated pathology, fully automated segmentation algorithms for multiple 
layers may be too time-consuming for routine use and manual user intervention 
may not be practical. For example, in the case of an RPE detachment, the layers 

1 5 that are normally brightly reflecting very near to the reference surface are now 

very far removed from the reference surface. The human eye is very sensitive to 
hue changes, and a convenient way to illustrate such RPE detachment is to use 
hue changes. Fig. 4 shows a 2-D retina image with the hue mapped as the square 
of the distance from the fitted RPE reference surface. 

20 (4) a flattened 3-D rendering of a retinal surface or surfaces, which shows the 

elevation relative to a reference surface rather than its actual contour elevation in 
the image data can provide a more meaningful view of the retina anatomic 
features. Fig. 7 shows a three-dimensional rendering of the ILM elevation relative 
to a paraboloid fitted RPE reference surface. Different pupil positions cause tilt in 

25 the recorded retinal images, and variations in working distance cause different 

curvature in the recorded retinal images. Warping the OCT data to flatten this 
image may have advantages for more standardized presentation, regardless of 
exact pupil and z position, aiding comparisons of the images between visits, or 
registration of multiple scans for other purposes such as speckle reduction. 

30 (5) a 3-D rendering of the ILM with a color mapping for elevation relative to the 

reference surface fitted to the RPE can combine the false color image previously 
described in (1), along with the actual retina thickness which could indicate the 
presence of traction by membranes on its surface. Fig. 8 shows a three- 
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dimensional rendering of the actual retina thickness superimposed with a 
pseudocolor image indicating ILM elevation with respect to the paraboloid fitted 
RPE reference surface. On the other hand, such a superimposed 3D rendering can 
also be the other way round. For example, Fig. 9 shows a three-dimensional 
rendering of the ILM surface elevation relative to the paraboloid fitted RPE 
reference surface superimposed with a pseudocolor image indicating actual retinal 
thickness. The ILM elevation may reflect the position in the image data or a 
rendering that is flattened to a reference surface as previously described in (4). 
[0036] Axial resolution may be wasted if the z-range of the scan does not follow 

the contour of the retinal tissue. A few initial scans could be used to determine the 
reference surface, then a retina-following scan could be performed by changing the OCT 
reference arm length to follow the predetermined reference surface as the transverse scans 
are performed. 

[0037] Note that the present invention can also be applied to B-scan images in 

which case, the term reference surface should interpreted as a reference curved line and 
tissue layer/boundary will also be interpreted as a curved line. Fig. 4 shows a B-scan 
retina image with the hue mapped as the square of the distance from the fitted RPE 
reference surface. The vertical distance in the image relative to the location of the fitted 
RPE reference surface is encoded as a color which is used to highlight the image. 
[0038] The present invention does not need to follow the exact sequence as shown 

in Fig. 3, as other additional steps can be inserted to perform substantially equivalent 
operations. For example, the fitting operation may be approximated by smoothing or 
otherwise filtering the retinal layer. Also, the reference surface might not be the direct 
result of fitting, but some filtered version thereof Furthermore, a variation of this idea 
could use two such reference surfaces, rather than the elevation from an unfitted surface 
to a reference surface. 

[0039] The presently invented method could be applied to the analysis of the 

retina or curvature of the eye in existing and future OCT systems. It can also be used for 
analysis of other biological tissues such as the skin. Also, it may find use in ultrasound 
and confocal microscopy systems as well. 

[0040] Fig. 1 1 shows a basic spectrometer based spectral domain OCT system 

1 100. The light wave from the broadband emitter 1 1 10 is preferably coupled through a 
short length of an optical fiber 1 1 12 to an input port (port I) of a fiber optic coupler 1114, 
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which splits the incoming light beam into the two arms of a Michelson interferometer. 
The two arms each have a section of optical fiber (1116 and 1118) that guides the split 
light beam from the two output ports (port II and port III) of the fiber coupler 1 1 1 4 to a 
sample 1 124 and a reference reflector 1 126 respectively. For both the sample arm and the 
5 reference arm, at the terminating portion of each fiber, there may be a module containing 
optical elements to collimate or focus or scan the beam. Illustrated in Figure 1 1 as an 
embodiment are two focusing lenses 1 120 and 1 122. The returned light waves from the 
sample 1 124 and the reference reflector 1 126 are directed back through the same optical 
path of the sample and reference arms and are combined in fiber coupler 1114. A portion 

10 of the combined light beam is directed through a section of optical fiber 1130 from port 
IV of the fiber coupler 1 1 14 to a spectrometer 1 1 50. Inside the spectrometer, the light 
beam is dispersed by a grating 1 152 and focused onto a detector array 1 154. Note that the 
principle of operation of a tunable laser based swept source OCT is very similar to that of 
a spectrometer based spectral domain OCT system (see for example, Choma, M. A. et al. 

1 5 (2003). "Sensitivity advantage of swept source and Fourier domain optical coherence 

tomography." Optics Ex press 11(18): 2183-2189), hence the spectral domain OCT system 
for obtaining the 3D image data set can also be a swept source OCT system. 
[0041] Although various embodiments that incorporate the teachings of the 

present invention have been shown and described in detail herein, those skilled in the art 

20 can readily devise many other varied embodiments that still incorporate these teachings. 
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1 . A method for generating images from a 3D image data set, said data set 
obtained from scanning a sample, said sample having at least one boundary surface 
associated therewith, said method comprising the steps of: 

identifying a subset of the image data corresponding to the boundary 
surface; 

fitting the subset of image data to a substantially smooth reference surface; 
calculating the distance between points on the reference surface to points 
of interest within the sample identified from the 3D image data set; and 
generating an image map of the calculated distances. 

2. A method as recited in claim 1 , wherein the points of interest correspond 
to points on said boundary surface. 

3. A method as recited in claim 1, wherein the points of interest correspond 
to points on a different boundary surface. 

4. A method as recited in claim 1 , wherein the sample is skin. 

5. A method as recited in claim 1 , wherein the sample is a blood vessel. 

6. A method as recited in claim 1, wherein the sample is an eye. 

7. A method as recited in claim 6, wherein the boundary surface is the 

cornea. 

8. A method as recited in claim 6, wherein the boundary surface is the retina. 

9. A method as recited in claim 1, wherein the 3D image data set is obtained 
using an OCT system. 
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1 0. A method as recited in claim 1 , wherein the 3D image data set is obtained 
from a confocal microscope. 

11. A method as recited in claim 1 , wherein the 3D image data set is obtained 
from an ultrasound scanner. 

1 2. A method for generating information from a 3D image data set, said 3D 
data set obtained from scanning a sample, said sample having at least one boundary 
surface associated therewith, said method comprising the steps of: 

identifying a subset of the image data corresponding to the boundary 
surface; 

fitting the subset of image data to a substantially smooth reference surface; 
calculating the distance between points on the reference surface to points 
of interest within the sample identified from the 3D image data set; and 
generating a 2D data set of the calculated distances. 

13. A method as recited in claim 12, further including the step of generating an 
image map from the 2D data set. 

14. A method as recited in claim 12, wherein the 3D image data set is obtained 
using an OCT system. 

15. A method as recited in claim 14, wherein the sample is an eye. 

1 6. A method as recited in claim 15, wherein the points of interest correspond 
to points on said boundary surface. 

1 7. A method as recited in claim 1 5, wherein the points of interest correspond 
to points on a different boundary surface. 

1 8. A method as recited in claim 1 5, wherein the boundary surface is the 

cornea. 
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19. A method as recited in claim 15, wherein the boundary surface is the 

retina. 

20. A method as recited in claim 12, wherein regions of the boundary surface 
that are significantly impacted by disease are omitted from the fitted surface to minimize 
perturbations of the fitted surface. 

21. A method for generating images from a 3D image data set, said data set 
obtained from scanning a sample, said sample having at least one boundary surface 
associated therewith, said method comprising the steps of: 

identifying a subset of the image data corresponding to a first boundary 
surface; 

fitting the subset of image data to a substantially smooth reference surface; 
identifying a second boundary surface; 

calculating the distance between points on the reference surface to 
corresponding points on the second boundary surface; and 
generating an image map of the calculated distances. 

22. A method as recited in claim 21 , wherein the first and second boundary 
surfaces are the same. 

23. A method as recited in claim 2 1 , wherein the first and second boundary 
surfaces are different. 

24. A method as recited in claim 2 1 , wherein the 3D image data set is obtained 
using an OCT system. 

25. A method as recited in claim 24, wherein the sample is an eye. 



26. 

cornea. 



A method as recited in claim 25, wherein the boundary surface is the 
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27. A method as recited in claim 25, wherein the boundary surface is the 

retina. 

28. A method as recited in claim 21 , wherein regions of the first boundary 
surface that are significantly impacted by disease are omitted from the fitted surface to 
minimize perturbations of the fitted surface. 

29. A method for generating an elevational map from a 3D image data set, said 
data set derived from an OCT scan of an eye, said eye having at least one boundary 
surface associated therewith, said method comprising the steps of: 

identifying a subset of the image data corresponding to the boundary 
surface; 

fitting the subset of image data to a substantially smooth reference surface; 
calculating the elevation between points on the reference surface to points 
of interest within the sample identified from the 3D data set; and 

generating an elevational image map from the calculated elevational 

points. 

30. A method as recited in claim 29, wherein the points of interest correspond 
to points on said boundary surface. 

31. A method as recited in claim 29, wherein the points of interest correspond 
to points on a different boundary surface. 

32. A method as recited in claim 29, wherein regions of the boundary surface 
that are significantly impacted by disease are omitted from the fitted surface to minimize . 
perturbations of the fitted surface. 

33. A method as recited in claim 29, wherein the fitting step is performed by 
one of a second-order polynomial fit, Zernike polynomials, Chebyshev polynomials, 
Bessel functions. 
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34. A method as recited in claim 29, wherein the subset of image data is fit to 
portion of a sphere. 

35. A method as recited in claim 29, wherein the fitting step is performed 
using a low-order polynomial fit in x and y. 

36. A method as recited in 29, claim wherein the reference surface is created 
by fitting to 3D data corresponding to the retinal pigment epithelium. 

37. A method as recited in claim 29, wherein the elevational points are 
calculated between the reference surface and the retinal pigment epithelium. 

38. A method as recited in claim 29, wherein the elevational points are 
calculated between the reference surface and the inner limiting membrane. 

39. A method as recited in claim 29, wherein the elevational information is 
presented as a contour map. 

40. A method as recited in claim 29, wherein the elevational information is 
presented as a color map. 

41 . A method as recited in claim 29, wherein the elevational information is 
presented as 3D surface rendering. 

42. A method as recited in claim 29, further including the step of generating 
secondary image information related to a structural feature of the eye and wherein said 
map comprises both the elevational points and the secondary image information. 

43. A method as recited in claim 42, wherein the elevational points are 
presented in a form selected from the group consisting of contours, color variations, grey 
scale, 3D surface rendering and wherein the secondary image information is presented in 
a form selected from the same group as the elevational points but using a different form 
than used for the elevational points. 
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44. A method as recited in claim 42, wherein the elevational points are 
presented as a color map and the secondary image information is shown as a 3D surface 
rendering. 

45. A method as recited in claim 42, wherein the elevational points are 
presented as 3D surface rendering and the secondary image information is presented as a 
color map. 

46. A method as recited in claim 42, wherein the secondary image information 
corresponds to one of the elevation between actual boundaries, actual thickness of a tissue 
layer, image signal strength, birefringence and data from a fit to a different reference 
surface 

47. A method as recited in claim 42, wherein the elevational points are 
presented as a color map and said secondary image information is image signal strength 
and is presented as variations in brightness in the color image 

48. A method as recited in claim 42, wherein the image map comprises a 
contour map having variations in color superimposed thereon with the elevation points 
corresponding to either the contour or color variations and with the secondary image 
information corresponding to remaining image variable. 

49. A method as recited in claim 42, wherein the boundary surface is the RJPE 
and wherein the secondary information corresponds to the distance from the ILM to the 
RJPE and is displayed in the form of a contour map and wherein the elevational points 
correspond to the distance from ILM to the reference surface and are displayed in the 
form of a color map superimposed on the contour map 

50. A method for generating an elevational map from a 3D image data set, said 
data set derived from an OCT scan of an eye, said eye having at least one boundary 
surface associated therewith, said method comprising the steps of: 



WO 2007/028531 



18 



PCT/EP2006/008465 



identifying a subset of the image data corresponding to the boundary 
surface; 

fitting the subset of image data to a substantially smooth reference surface; 
calculating the elevation between points on the reference surface to points 
of interest within the sample identified from the 3D data set; and 
generating a 2D set of the calculated elevational points. 

51. A method for generating images from a 2D image data set, said data set 
obtained from scanning a sample, said sample having at least one boundary surface 
associated therewith, said method comprising the steps of: 

identifying a subset of the image data corresponding to the boundary 
surface; 

fitting the subset of image data to a substantially curved reference line; 
calculating the distance between points on the reference line to points of 
interest within the sample identified from the 2D image data set; and 
generating an image map of the calculated distances. 

52. A method as recited in claim 51 , wherein the image data set includes 3D 
information. 

53. A method for generating images from a 2D image data set, said data set 
obtained from scanning a sample, said sample having at least one boundary surface 
associated therewith, said method comprising the steps of: 

identifying a subset of the image data corresponding to the boundary 
surface; 

fitting the subset of image data to a substantially curved reference line; 
calculating the distance between points on the reference line to points of 
interest within the sample identified from the 2D image data set; and 
generating a 2D data set of the calculated distances. 

54. A method as recited in claim 53, wherein the image data set includes 3D 
information. 
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55. A method as recited in claim 54, further including the step of generating an 
image map from the 2D data set. 
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